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We present rigorous bounds on the thermalization time of the family of quantum mechanical spin systems 
known as stabilizer Hamiltonians. The thermalizing dynamics are modeled by a Davies master equation that 
arises from a weak local coupling of the system to a large thermal bath. Two temperature regimes are considered. 
First we clarify how in the low temperature regime, the thermalization time is governed by a generalization 
of the energy barrier between orthogonal ground states. When no energy barrier is present the Hamiltonian 
thermalizes in a time that is at most quadratic in the system size. Secondly, we show that above a universal 
critical temperature, every stabilizer Hamiltonian relaxes to its unique thermal state in a time which scales at 
most linearly in the size of the system. We provide an explicit lower bound on the critical temperature. Finally, 
we discuss the implications of these result for the problem of self-correcting quantum memories with stabilizer 
Hamiltonians. 


The study of thermalization times of Markov processes for 
classical spin systems is by now a mature research field, which 
has lead to important contributions for stochastic optimiza¬ 
tion problems, image reconstruction and a better understand¬ 
ing of Monte Carlo algorithms. Aside from these applications, 
the analysis is interesting in its own right as it has lead to a 
better understanding of non-equilibrium thermalization pro¬ 
cesses typically modeled by Glauber dynamics ll]]. In this let¬ 
ter, we study quantum mechanical Markov processes that con¬ 
verge to the Gibbs state of interacting quantum spin Hamil¬ 
tonians. We are in particular interested in the equilibration 
time of a specific type of master equation, the Davies gener¬ 
ator ||2l, which can be seen as a physically motivated quan¬ 
tum generalization of classical Glauber dynamics. We prove 
thermalization time bounds for the important class of stabi¬ 
lizer Hamiltonians. These models are well motivated by the 
study of self correcting quantum memories (SCQM) and the 
bounds provide rigorous no go results for certain candidate 
models. We show that at low temperatures the convergence is 
determined by a generalized energy barrier (Eqn. (HI) of the 
stabilizer Hamiltonian - a quantity that only depends on the 
excitations of the Hamiltonian. Furthermore, we give a bound 
on the critical temperature above which the thermal process 
always equilibrates rapidly and is independent of the low en¬ 
ergy excitations. 

We bound the thermalization time of a quantum dynami¬ 
cal semigroup generatored by a Liouvillian C in terms of the 
mixing time. If the steady state p of £ is the Gibbs state of 
a Hamiltonain H, then the mixing time tmix is the least time 
such that \ \ip{f) — p||i < e“^, for any initial state tp and t < 
tmix, where ip{t) = (<p), and the semigroup converges 

to the Gibbs state p = e“^^/tr Here and in the 

following, C* refers to the dynamics in the Schrodingier pic¬ 
ture, and to C in the Heisenberg picture (see Eqn. (|2]i). It can 
be shown in general ([3l, that ||:p(f) — p||i < 
where ||p“^|| < is the inverse of the smallest eigen¬ 

value of p and A is the spectral gap of the Liouvillian C, de¬ 
fined as the smallest (in magnitude) non-zero real part of an 
eigenvalue of C. Since ||iT|| scales as the volume of the sys¬ 


tem, the mixing time primarily depends on the gap of the Li¬ 
ouvillian: tmix ~ |iT| I). Hence our main theorems 

will consist of lower bounding the spectral gap of the genera¬ 
tors C in different settings. 

A central motivation for this work is the investigation of 
thermal stability of candidate Hamiltonians that can serve as a 
SCQM 101. Many of the quantum spin Hamiltonians we con¬ 
sider are topologically ordered ||3l and can encode quantum 
information in a degenerate subspace, for which the states can 
not be distinguished by local operations. Ideally, a SCQM 
should be able to store an arbitrary quantum state in contact 
with a thermal bath that survives for a time which grows ex¬ 
ponentially in the size of the system below some critical tem¬ 
perature. In other words, the code subspace would need to 
be metastable. So far, only one candidate system is known 
which satisfies this definition: the 4D toric code 0 . In two 
or three dimensions a number of no-go results exist IMS, 
which exclude the possibility of a quantum memory for cer¬ 
tain classes of Hamiltonians. Most of these results rely on 
heuristic assumptions, such as the validity of the Arrhenius 
law d . The bounds on the thermalization time in this paper 
provide rigorous no go results, since tmix provides a natural 
upper bound to the lifetime. If the system thermalizes in a 
time that scales polynomially in the system size, we will say 
the system is not a good quantum memory. Note that this 
approach only permits to rule out possible candidates since a 
slow thermalization rate can only indicate the ability to store 
classical information, while quantum information may have 
already been lost. 


Preliminaries: The stabilizer formalism was introduced by 
Gottesman 11211 as a convenient framework to express quan¬ 
tum error correcting codes (QECC). It soon became clear that 
this formalism could also express a large family of commut¬ 
ing, frustration-free Hamiltonians, with the toric code as a 
prime non-trivial example fllH]. The Pauli group on N qubits 


is defined as the set Vn = ® ... ® off' : ay € 

{I, X, Y, Z}, ip G {0, 7r/2, TT, 37r/2}} with the group oper¬ 
ation being the usual matrix multiplication and where {cr“^} 
are the Pauli operators at site j. The group Vjy has 4 ^+^ 
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elements, half of these are Hermitian (and the other half anti- 
Hermitian). A Pauli stabilizer group S is an Abelian sub¬ 
group of Vn such that —I ^ S. By excluding —I from S, 
we guarantee that all the elements of s G 5 are Hermitian 
and can stabilize a subspace with eigenvalue (+1). The sub¬ 
space C of vectors |'0) satisfying Vs G S, s|^) = |'0) is called 
the subspace stabilized by the group S. To define the stabi¬ 
lizer Hamiltonian associated to S we need to introduce the 
generating set Q C S of the stabilizer group. Recall that ev¬ 
ery Pauli stabilizer group may be written as the group gen¬ 
erated hy M < N hermitian and commuting Pauli operators 
Q = {gi,..., Pm}- If the generators are all independent the 
corresponding stabilized subspace C is dimensional. 

A stabilizer Hamiltonian iV is a Hermitian operator on (D^ 
which may be written with constants Jk > 0 as 

M 

H = -J2Jk9k ( 1 ) 

k=l 

of Stabilizer operators gk G G- Note that the Hamiltonian 
is only defined in terms of the generating set G and not the 
group S. The excitations therefore depend on the choice 
of generators. The ground space, can be identified with the 
stabilized subspace C. Stabilizer Hamiltonians are frustration 
free and can be diagonalized easily because all gi~ commute. 
Nevertheless they can exhibit exotic quantum behavior, such 
as topological order. For a review and a class of example 
Hamiltonians, the reader is referred to Ref. 0. 

The Davies generator: We model the thermalizing dynam¬ 
ics of the stabilizer Hamiltonians by weak (or singular) lo¬ 
cal couplings to a large memoryless bath in thermal equilib¬ 
rium. In this setting, there is a standard procedure for de¬ 
riving a Markovian master equation describing the dissipative 
action of the bath on the system 101. Without loss of gener¬ 
ality, we can assume that the system-bath interaction Hamil¬ 
tonian is given by = J2k,ak ® where 

are some set of operators acting on the bath and Sa^, = 
are the system coupling operators acting on a single site as 
Pauli operators, which span the full local matrix algebra. 
The master equation that is obtained from this procedure is 
£(/) = i[Hes, f] + D{f), where / is an observable, and Hgg 
is an effective Hamiltonian which does not contribute to the 
spectrum of C. The map D describes the dissipative part of 
the evolution, and is given by 


nf) = E , (2) 

ockt^ ^ ' 

where ui = Ei — Ej are the Bohr frequencies of H with eigen¬ 
values Ei. The rate with ha^{—uj) = e~^^ha^{uj) 

is determined by the bath auto-correlation function 
and encodes the dependence on /3. We assume that 
hmin = mine is the smallest transition rate over 

all realized Bohr frequencies. The jump operators are the 


Fourier coefficients of the time evolved system coupling op¬ 
erators exp{iHt)Sak exp(—iiJf) = - 

{u!,ak) the action of the jump operators can be understood as 
mapping eigenvectors of H with energy E to eigenvectors of 
H with energy E + u!, and transfer energy ui from the system 
to the bath and back. The generator satisfies quantum detailed 
balance with respect to the Gibbs state p of the Hamiltonian 
H 101. This ensures that p is a fixed point and furthermore 
that the spectrum of £ is real. 


The low temperature regime: It has been shown 11 1311 that 
Davies generators of bounded one dimensional lattice systems 
are always gapped. The proof hinges on an equivalence the¬ 
orem between the relaxation time of Gibbs samplers and the 
correlation properties of the Gibbs state. Hence, in ID there 
exists a correspondence between the static and the dynamical 
critical behavior. We want to identify properties of the Hamil¬ 
tonian that lead to a short mixing time independent of temper¬ 
ature, in dimensions larger than one. The following can be 
shown: the central quantity which determines the spectral gap 
at low temperature is a slight generalization of the well known 
energy barrier for logical operators | 0 ]. 

In order to define the generalized energy barrier e, we need 
to choose an enumeration F of the lattice points, by which 
every Pauli operator is constructed by single qubit operations. 
For any p G 7 ^at, we construct a path in Pauli space by setting 
po = 1 and defining at every step I G F the Pauli operator gi 
to correspond to 77 up to site t and identity on the rest, so that 
the full Pauli is 77 = 77 ;^ with the largest element in F. We 
define the reduced generating set of the Pauli operator 77 by 
Gti = {g & G [g,v] = 0}, i.e. all terms in the Hamiltonian 
that commute with 77 . If we set J = max^ \Jk\ the energy 
penalty of 77 w.r.t. to the ordering F is defined as 


er(77) = 2Jmax#{5 G ; [gug]^'^)- (3) 

The constant that determines the spectral gap is now obtained 
by considering the highest energy penalty of any Pauli string 
with respect to the best labeling of lattice sites. That is we 
define 


e = min max erf??)- 
r ^dVN 


(4) 


This quantity gives rise to the following spectral gap bound of 
the Davies generator. 


\d > 


mm ^-2l3e 


4£ 


(5) 


where £ = |r| denotes the length of the longest Pauli path 
This result was first shown in 11411 . A bound to Ad is ob¬ 
tained by a direct evaluation of the Poincare inequality (Supp. 
Mat. Eqn. (l28T l) for the Davies generator. The gap can be 
related to the energy barrier using methods developed in Refs. 

HS. 

The bound in Eqn. (|5]l is stated in terms of the optimal 
choice of F. However, any non optimal ordering also gives 
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FIG. 1. (color online) For illustration purposes, we focus on the two 
simplest example systems, the toric code (TC) fig. (a), to explain the 
low temperature bound and the Ising model (IM) fig. (b) for the high 
temperature bound. Figure (a): The TC (1 describes a lattice system 
of N spins at every link on the lattice indicated by black / white 
dots. The Hamiltonian H — — J — J Bp is given as a 

sum of (blue) vertex terms Ay = cr^ cr^ for every vertex v 

and (grey) plaquette terms Bp = cr^ cr^ for every plaquette 

p, so that Q — {Ay, Bp}y^p. We have depicted an example Pauli 
p G Vn by (red) X, Y, Z’s placed on the links. The corresponding 
reduced set Qp has also been depicted by removing the appropriate 
(blue) vertex and (grey) plaquette terms from Q. Figure (b): The IM 
is understood as a stabilizer Hamiltonian where the qubits reside on 
the vertices : H = —erf cr^,so that Q = {af We 

have depicted the set of qubits in Nj for some j as a grey rhombus 
and chose a clockwise order of the qubits. We have IWjl = 5 and 
Sj = {ctf ctf}ieNj\{j} SO that \Sj \ = 4 and furthermore IS'inSjl = 
1. We can evaluate Eqn. ® as k{I3) = — 1) from 

which ~ 249. 


rise to a valid bound. The quality of the bound strongly de¬ 
pends on the choice for T. For the 2D toric code, for ex¬ 
ample, a good choice gives rise to e < 2 J, whereas a poor 
choice can easily lead to an exponentially small bound on An 
with e < 0{y/N). For a suitable ordering of the lattice sites 
F, one typically has that the longest path = 0{N). We 
expect that the dependence on in Eqn. (|5ll is an artifact 
of the derivation, and that the true bound should rather be 
Xd > In particular in the limit /3 —>■ 0 Eqn. 

dSll predicts A = whereas the true value of the gap 

is known to be A = 0(1) ifisll . If the Hamiltonian is not 
protected by an extensive generalized energy barrier, i.e. we 
have e = const, the system thermalizes in a time that scales 
as tmix = 0{N‘^) since ||iT|| = 0{N) for local stabilizer 
Hamiltonians. 

We now consider a simple example of er(?7) for the 
situation illustrated in Eig. [Ha). We choose an ordering F 
where all lattice points are traversed twice as to first build up 
the Z and then the X factors of rj. The Y Pauli’s are read as 
first applying Z and then X on the same site. We have that 

= |r| = 2N. The qubits are ordered so that for Z factors 
we start to traverse the white qubits in each column (top to 
bottom) moving to the next column on the right. Then we 
traverse each row of black dots (left to right) and move to the 
next lower row. For the W-factors the roles are reversed. We 


first order the black qubits in each column (top to bottom) 
and then traverse the rows of white qubits (left to right). 
The ordering fixes the path pt for any Pauli matrix p G Vn- 
This amounts to decomposing every p into small patches of 
incomplete ’’logical operators” |@1. We see that along this 
path at most one G Grj is violated so that er{p) < 2 J, for 
all p. This implies a lower bound Xd > 

Had we reversed the ordering F for X and Z the bound would 
be exponentially small. The generalized energy barrier e is 
related to the standard energy barrier of logical operators 
101 when the lattice can be decomposed into the support of 
equivalent minimal logical operators. Then the argument 
above can be generalized and e corresponds to the largest 
energy barrier of the canonical logical operators. 


The high temperature regime: It is expected that the low 
energy excitations do not determine the equilibration times at 
sufficiently high temperature, since this behavior is typically 
associated with the disordered phase of the model. In Ref. 
lll3ll . it was shown that there must exist a transition tempera¬ 
ture above which the spectral gap of the Liouvillian is constant 
in the system size. From physical intuition one would expect 
that the transition temperature that gives rise to a constant gap 
in the dynamics should coincide with the static critical temper¬ 
ature of the Gibbs state of the spin Hamiltonian iflbl] . That this 
is the case has only been shown in specific cases Cl [Hill 
and it is not known whether it holds for general quantum spin 
models. We therefore only discuss the dynamical critical tem¬ 
perature for which rapid mixing sets in here. Below we give 
an explicit lower bound on this temperature. 

The bounds only depend on the coordination properties of 
the Hamiltonian and can be evaluated by considering the local 
neighborhood of single spins. The neighborhood Nj of spin 
j is defined as the set of qubits I G Nj that share a stabilizer 
generator, so that l,j G supp(g) for some g G G- This set 
essentially coincides with the support of all the . We are 
free to choose an ordering of the qubits in Nj so that j = 1 is 
the smallest element in this order. Furthermore we define the 
set Si C G SiS, the set of generators that have support on site /, 
i.e. g G Si if I G supp(p). In Proposition (Sup. Mat. Prop. |3 
we prove a lower bound on a constant n that can be simplified 
to 


k(/?) = max3(|iV,-| - ^ 

^ iGNj k>l 

( 6 ) 

Here J = maxfcg^^^ Jk is the largest coupling constant. 
The bounds on temperature and gap are obtained from the fol¬ 
lowing theorem. 


Theorem 1 Above a threshold temperature /?*, the gap of C 
is independent of the system size. We get for any /I < /3, := 
max{/l|K(/3) < 1}, that Xd > (1 — n{P)), 

with S', = miuj |Sj|. 


Proof. 

i3. 


The proof follows closely ideas developed in Ref. 
We defer the technical calculations of constants to 
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the appendix, and only present the main steps here. To 
start, we define the so-called heat-bath generator: Q*{cr) = 
Ej = E, - a, where E*(a) = 7 ]tr, (a) 7 ,, 

and 7 = (trj here tr^ (/) denotes the 

normalized partial trace and corresponds to a single qubit 
twirl. Q is a proper Liouvillian, and by a comparison the¬ 
orem (Supp. Mat. Lemma|3ll we can show that the gap of the 
Davies generator can be lower bounded by the gap Xq of Q 
as: 


Xd > RXq, (7) 

where we have the bound R > which is system 

size independent. Hence, it suffices to lower bound the gap of 
the heat-bath generator, to get a lower bound on the gap of the 
Davies generator. The spectral gap of Q can be bounded by a 
generalization of the Aizenman and Holley d conditions to 
quantum spin systems H. 

Let Tj = e*® be the heat-bath semigroup in the Heisenberg 
picture (i.e. acting on observables), and define the discrete 
gradient: dkif) = f - ti'fc (/). Note that dkQkif) = dk{f), 
since dk^kif) = 0. Therefore, 


-^dkTsif) = -dkUf) + J2Qjdk{f) + J2[dk,Qj]Uf) 


If we define = exp(t Ej^k )’ 


ds 


(e^T.^lldkTsif)) = Y, e^Tt_s[dk, Q,]T,(/) (8) 


Now, integrating from 0 to f on both sides of Eqn. (|8]l 
and using the fact that the semigroup T* is contractive (i.e. 
||'rt(/)|| < ll/ll, for all /), we get 


\\dkTt{f)\\<e-*\\dkf\\ 


E 




\\[dk,QMm 


The trick now is to bound the last term point wise using 
Lemma (Supp. Mat. Lemma|4]i, which claims that there exist 
a set of positive constants a* > 0 such that 

\\[dk,Qm\\<Y<k\\di{f)\l (9) 

I 


with J2k Ej^k fc — where k < 00 for any I, k. Then we 
get: 

|im(/)|||<e-‘|||/||| + « /‘e-(‘-^)|||r,(/)|||, (10) 

Jo 

where lll flljl = Ek ll^fe/ll ’^^e so-called oscillator norm 
ifisi 1^ 1^ . Iterating this equation yieldls H|r^(f)lll < 
g-(i-«)‘|||j||| A telescopic sum argument fllSll can then be 
applied to obtain 

im(/)-tr(p/)||<e-(i-'=)‘|||/|||, (11) 


where p is the Gibbs state of H. Clearly, whenever k < 1 
we get a non-trivial bound. Invoking properties of Lp norms 
(Supp. Mat. Lemma |6l), Eqn. (fTTI) implies that the spectral 
gap of Q is lower bounded by Aq > 1 — k. The final step in 
the proof consists in finding an upper bound on k. We obtain 
the bound Eqn. (|6]l from (Sup. Mat. Prop. |5ll after additional 
simplifications, so that applying Eqn. (|7]i yields the result. □ 


Typically we will find that <C 1 so that we can consider 
the expansion of k and only focus on the first order term to 
obtain a rough estimate of /3,. We have that Ri 6^ J (| TVj | — 
1) EieN- Ek>i 1‘^i ^ ^hat we can immediately solve 

At = 1 for /3,. 

The bound we obtain yields the same functional depen¬ 
dence on the temperature as in the static case 0. However, 
it is clearly not optimal in practice. Eor instance, for the 
Ising model in two dimensions, our bound yields a dynamical 
critical temperature of J)“^ ^ 249, whereas the correct 
static critical temperature is ~ 2.27. Eor comparison, in Ref. 
IB a bound on the static critical temperature corresponded 
to a value of ^ 24.58. Since the bound only depends on 
the local properties of the stabilizer Hamiltonian it can be 
evaluated easily, however this implies that the long range 
aspects, needed for a tight bound on the critical temperature, 
are neglected. The fact that for /? < /3, we have that 
A = const implies again with ||iJ|| = 0{N) for local 
stabilizer Hamiltonians that tmix = 0(N). 


Conclusions: We have presented a complete characteriza¬ 
tion of the thermalization time for stabilizer Hamiltonians. 
By now it is known iB 23] that the existence of an energy 
barrier is not sufficient for the thermal stability of a SCQM. 
However, the low temperature result shows that at least for 
Pauli stabilizer Hamiltonians the existence of an energy 
barrier e is a necessary condition. This implies in particular, 
that mechanism such as the ones analyzed in 11241251] . where 
random perturbations in the couplings suppress the coherent 
propagation of excitations, cannot protect against thermal 
errors. Eurthermore, we find an explicit bound on the critical 
temperature above which the gap of the generator is constant 
in the system size. The bound confirms the intuition that 
at sufficiently high temperature the particular low energy 
properties of the model do not influence the thermalization 
dynamics and only local properties of the model matter. 
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SUPPLEMENTAL MATERIAL 
Stabilizer Heat-bath and Davies generators 

In this subsection, we will analyze the Heat-bath Liouvillian Q defined in the proof of Theorem[T]and derive a more explicit 
representation for stabilizer Hamiltonians. Moreover, we will introduce notation that will facilitate the derivation of the constants 
discussed in the main body of the text. Since the terms in the Hamiltonian 

H = -Y,Jkgk ( 12 ) 

k 

commute, it can be diagonalized in the common eigenbasis of all {gk}- The individual gk have eigenvalue ±1, so the projectors 
onto their eigen-subspaces are Ilk{bk) ■= 5(1 + (—l)^''pfc). Here the choice bk = 0,1 denotes the projector onto the local 
subspace with eigenvalue (—1)^'“. Since all {gk} commute, we can construct global projectors by taking their products. Hence, 
eigenvalues and corresponding subspaces can be described by an M = \Q\ component bit string b. That is, we have that for 
every b G {0,1}^, the eigenvalues are 

e{h) with projectors P{h) := Ui{bi)U 2 {h 2 )...UM{bM)- (13) 

k 

With this at hand we can diagonalize the Hamiltonian immediately and obtain 

H = Y,^ib)P{b). (14) 

b 

The reader familiar with the stabilizer formalism will immediately recognize the bit string b as the syndrome of the stabilizer 
code d generated by Q. The Gibbs state of H can immediately be written as Z = Z ^'^^^'^P{b). For 

convenience we will also frequently write pa = Z~^ exp(—/?e(a)), so that p = Z~^e~^^ = PaP{ot)- 

To analyze the action of the Davies or Heat-Bath generators, it is necessary to understand how the Pauli operators Vff act on 
the eigen projectors of the Hamiltonian. We therefore need to introduce some notation. Recall that the Pauli group is isomorphic 
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to Vn — 22^“''^ = {0,where two bits are needed to encode the phase information. Since we work with the Pauli 
matrices as a basis we do not need to carry the phase information and can identify every Pauli operator up to a phase with a 2N 
- bit string a € We write for this Pauli matrix a {a) € Vn- In particular we denote the single qubit Pauli operators acting 
on site j by a{aj) = (t“, where aj = (0,0)^, (0,1)^, (1, 0)j, (1,1)^ corresponds to 1 , cr^ respectively. The full algebra 

is indexed by concatenating these bit strings from every site. Given the set Q of Pauli’s which comprise the Hamiltonian H, we 
can now define a syndrome, or excitation to every Pauli cr(a) to which we will refer as e(a) € TlJ^ which is a M-bit string 
defined as e(a) = (ei(a, )... eNrict)) defined by commutation properties a{a)gk = (—for every k. Hence the 
action of a Pauli operator on a projector P{a) can be understood as 

a{a)P{a) = P{a(B e{a))a{a), (15) 

by direct observation. For notational purposes we will write the bit string a“ := a 0 e(a) so that P{a 0 e(a)) = P{a°‘) and 
e{a 0 e(a)) = e(a“) respectively. Furthermore, it is helpful to define the difference between two eigenvalues e(a) and e(a“), 
when one is obtained from the other by Pauli excitation a{a). That is we will frequently consider the Bohr frequency 


OJ 


‘(a) = e{a) - e(a“) = - 'V 2 Jfe(-l)“'=efc(a), 


fc 


(16) 


where by abuse of notation we read ek{a) G S 2 as a {0, l}-valued variable with standard addition in R. The purposes of 
the newly introduced notation is to encode the excitations of the stabilizers Hamiltonians in a compact manner that is model 
independent. This is approach was initiated in 11411 to obtain the low temperature bound in Eqn. (|5]l in the main body of the text. 
With this notation at hand, we are now in the position to analyze the action of the heat-bath and Davies maps of H on the matrix 
algebra <D[Vn] — M 2 n{<P). 


The heat-bath generator: We now analyze the heat-bath generator Q = Qj that was introduced for the proof of 
Theorem [T] The Heat-bath Liouvillian can be understood as a quantum extension of the classical Glauber dynamics. Recall 
that Qj{f) = IEj(/) — /, with ]Ej(/) = tr^- and jj = tr^ (e“^^) The conditional expectation E is a 

quantum channel (cot map). We will write the generator in Lindblad form as this will simplify its spectral analysis. For a qubit 
lattice system A, the heat-bath Liouvillian Q can be written in Lindblad form as 

2(/) E E 

jeA 


with jump operators = Yb \Gj{b)P{h)a{aj) where cr{a.j) are the single qubit Pauli matrices acting on site j and 


G,{h) = 



E 





-1/2 


(18) 


Observe that the normalized partial trace tr^ (•) over a single site j (local depolarizing map) can be expressed as a sum over 
Pauli matrices a{aj) = { 1 , crj^, tr]^, ctJ} as 


ti-j (/) 


“3 


(19) 


We will now turn to expressing Q in a specific basis of operators. First observe that. 


tG(e =JEE® 

aj b 


aj b 


( 20 ) 


Where the last equality follows from Eqn. (flSl l by adding e{aj) to all b and relabeling. We then get 
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7, =(tr, 

( \ - 1/2 

i E E j (^E 

= E(^Er”'''^'’“'^j P(6)j 
\ - 1/2 

I e-^f^^^^'>P{b) = ^Gj{h)P{h). (21) 



We have defined Gj (b) in (fTSl) and made use of the definition a;“i {b) = e(b) — e(&“i). Recall that the Liovilian is given by, 
Qjif) = tTj — / so we can write the trace as the twirl over the single site Pauli group as explained above. Hence, we 

may choose the Lindblad operators 


= ( 22 ) 

b 

Moreover, note that since for any operator / we have that trj (/) commutes with any other operator that is only supported on 
site j we obtain 

a cx. 

= tij- (e“^^) (e“^^) tr^ (e“^^) = 1. (23) 

This leads in the end to the decomposition of the generator Qj into Lindblad operators as given above. 

Davies Generator: For completeness we also briefly discuss the particular form of the Davies generators for the spectral 
decomposition of H introduced above. Recall that the Liouvillian is obtained by weakly coupling the system iJ to a thermal 
bath aB = Z~^ exp(—)> where Hb is the bath Hamiltonian. The system-bath interaction is given by Hint = ® 

where Sa are local Hermitian operators only acting on the system, that couple to the bath operators It can be shown 
13] that in the weak coupling limit the reduced evolution on the system is described by a Master equation of the form 
dtp{t) = —z[iTeff,p] + P*{p)- Here iTeff is the effective, lamb shifted, system Hamiltonian. This Hamiltonian does not 
contribute to the spectral gap of the generator ll2^ and we will therefore be ignore it in the following. 

We will now consider the Davies generators and express them for commuting Pauli Hamiltonians. Recall the definition of 
Davies generator as given in Eqn. (|2l) 


D(/)= .;.(/), 

ui,aj 


(24) 


where the individual terms are 


V, 


,;.(/) = KM i 


(25) 


We know that due to the KMS condition, the Fourier transform of the bath auto-correlation function satisfies (w) = 
(—(u). To state a bound that is independent of the specific Bath properties, we assume that the transition rates are 
bounded from below and above by a possibly temperature dependent constant independent of the system size; /imax > hf{uj) > 
Kin > 0. Moreover, recall that 



UJ 


( 26 ) 
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We choose the local perturbations on site j simply as all the single qubit Pauli operators 5“ = cr(aj). Since iJ is a stabilizer 
Hamiltonian, we can write its time evolution as exp{iHt) = exp(ie(6)t:)P(6) by means of the spectral decomposition in 

Eqn. (fT4li From which the individual 5'“(a;) can be read off directly 


:=J2S[uj-uj^^ib)]a{a,)Pib), 

b 


where he have defined the following 6 function 

= IJ 


for a; = 0 
else. 


(27) 


Comparing spectral gaps 


It is sometimes possible to bound the convergence behavior of one generator by that of another. In our case this has the 
following consequence. Suppose we obtain a bound on the spectral gap Ag of the Heat bath dynamics, as was done in the 
proof of Theorem [T] we can then use this bound to obtain a bound on the spectral gap of the Davies generator by means of a 
comparison technique. The general argument for such a comparison of spectral gaps can be made along the following lines. The 
Poincare inequality J^l gives rise to a variational characterization of the spectral gap of a Lindbladian C. 

AVarp(/) < f (/) V/gC^^, (28) 

where the two forms are defined as 

Varp(/) =tr(^p^/V^/) - |tr(p/)p and £(f) =-tr (^pipC{f )y (29) 


We refer to £{f) as the Dirichlet from and to Varp(/) as the variance. A direct evaluation of this inequality was carried out 
in lll4|] and yields the inequality on the gap as stated in Eqn. Q. For the high temperature bound we use this inequality as an 
intermediate step. If we can find a constant r such that £^{f) < t£^ {f), where £^{f) and £^{f) are the Dirichlet forms for 
Q and V respectively, we can infer a lower bound on the spectral gap Ad of the Davies generator in terms of Ag from a simple 
chain of inequalities. 

Lemma 2 Let Ag denote the spectral gap of the heat bath Liovillian Q (Eqn. f l77| )) associated to the commuting Pauli Hamil¬ 
tonian H = — spectral gap Ad of the corresponding Davies generator D (Eqn. ( 1251 )) can be bounded 

by 


Ad > 7?Ag. 


(30) 


The constant R is given by 


R> ^/inline (31) 

where in the lower bound we have defined /imin = min^c.^- h°‘(uj°‘^ (a)) and S'* = maxjgA |Sj |. 

Proof. The Poincare inequality, Eqn. (1281) . gives a variational characterization of the spectral gap. From the inequality we see 
that the gap is the largest constant A that satisfies AVarp(/) < £{f) for all / G (D^ . Note that both the heat bath generator Q as 
well as the Davies generator D have the same fixed point p = Z~^ exp(—/37f). This means that we can use the same variance 
for both inequalities. Hence by virtue of AgVarp(/) < £^{f) and applying lemma[2 we get 

AgVarp(/)<SQ(/)<rS^(/), V/G (32) 

so that have Ad > SAg with R = . We can directly bound t using: max^j- Gj{a) =4“^ maxo ^ ^ to 

obtain the stated bound on R. □ 

Since both maps are very similar in their support, we will be able to compare the Dirichlet forms of Davies Generator and 

the Heat bath dynamics locally. That is we will proceed to split the forms up in ever smaller subsystems and identify the worst 

bound among these simpler systems. This will then lead to a direct comparison between the two forms. 



9 


Lemma 3 Let £^{f) = —tr Q(/)^ denote the Dirichlet form of the heat-bath generator and let £^{f) = 

—tr ^p 5 P p^'D^f)^ denote the Dirichlet form of the Davies generator, for any f G (D^ .Then 

£Q{f) < T£^if), (33) 


where 


T = max 

j',a 



(34) 


Proof The fact that we are considering commuting Hamiltonians ensures, that both the Davies generator as well as the thermal 
course graining map can be written as a sum of local generators that are supported only in the neighborhood of a single site. We 
shall write in the following 


£Q{f) = Y,£f{f) as well as ^ (/), (35) 

jGA jGA 

with the local forms £^ {f) = —tr ^P^/V^2i(/)) and £f{f) = —tr p p^Vpf'^, and A is the full lattice. If we can 

find constants Tj so that the local terms are bounded as £'^ (/) < Tj£^{f) we know that the choice r = maxj Tj will provide 
the desired bound. To find bounds on the local Tj we follow an approach taken in Ref. and recast the problem into a semi 
definite inequality. We vectorize the / G (D^^, by the natural isomorphism that maps / —7> | /) G and write both forms 
in terms of hermitian Matrices £f{f) = (/ \£f\ /) and £p{f) = (/ \£p\ /)■ The local number Tj can now be bounded by the 
smallest number such that 


rpf - £f > 0 , 


(36) 


is positive semi-definite. We will see that i.f and £f are very similar in form, so we will be able to find good bounds on r, easily. 


We only derive the matrix £'^ because a similar derivation has already been performed in Ref. Ill4l] for the Davies generator. 
We will express the matrix entries of £f in terms the Pauli matrices normalized with respect to the Hilbert-Schmidt inner 


product. We have that the action of (( 7 ( 7 )) = —ypQj((T( 7 ))y/p defines the matrix We evaluate 


£fp{l)) = E P~pU{AfApaP)} - Af a(7)A“ j vY, 

OLj ^ ^ 

by using the representation of in Eqn. (l22l) . and by virtue of p = paP{o) we obtain 

'(o-(7)) = ^E E ^ 5 a,bP{l) {cj{aj)P{a)a{aj),ap)} P{m) 

-P{l)a{aj)P{a)<T{pP{b)a{apP{m)j 


otj a.b.L' 


(37) 


(38) 


If we now make use of the commutation relations Eqn. (fTSl) between cr(p) and P(a) we have that 

£?PiP) = iEE (Gf(a“0 + Pia)aP). (39) 

aj a ^ A 

Here we have defined the phase 0^.^^ = ±1 depending on wether the Pauli matrices commute or anti-commute 

a(aj)cr{j) = 6»«.,^cr(7)cr(aj)- 

In this form it is possible to read off the matrix £^ in the Pauli basis. As discussed previously, there is a correspondence 
between and Vm- We only care about the Pauli matrices and not the phase information of the group. We associate to 

every Pauli matrix, labeled by 7 G a vector in using the state | H) = P’ through 


I 7 ) = o-( 7 ) ® 1 |H). 


( 40 ) 
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Moreover, note that we can express the projectors as P{a) = 2~^ y ^„(—■ We can define a dual basis for 
each coset of the Stabilizer group S = {Q) in the full Pauli group Vn- Following 12611 we denote the coset Basis elements by 
I ) and define them by 


l«7o) = ^ E I 


(41) 




Here a corresponds to the syndrome of the generator set a, and 70 labels the representative of the coset. Following Ill4ll . we can 
express the matrix in this basis and read off the coefficients from Eqn. (l3^ to obtain 


^3^ = E E .70 («) I “70 ) (^70 I - .70 .70 I ^70 ) (^70 I' 

[7o]6Pn\S a^OLj 


(42) 


Here we have defined the following functions 


Al^{a) = - + Gf (a“'^)) ,/J^ and . 


(43) 


Note, that these functions are constant over the coset [ 70 ], i.e. one can easily verify that if we choose a different repre¬ 
sentative 7 i of the same coset we have that = A^^^{a) and = B^^^{a). Moreover the same holds for 

|o 7 o)(® 7 ol = l® 7 i)(® 7 il well as for 0 ( 13 , 70 1 070 ) (a“o I = (* 03 , 71 1 071 ) (““i I- This can be verified directly from the definition in 
Eqn. drill . Hence, the choice of representative is arbitrary. 


Essentially the same analysis for the Davies generator has already been done in Ref. Ill4|] to obtain the corresponding Dirichlet 
matrix. It turns out that the matrix is of the same form and we only need to replace A and B in Eqn. (l42l) by the corresponding 
versions for the Davies generator, which are 


'^a,7 


(a) 


\iK 


^a-ya-y ) s/Pa Pa-1 


and B^,^{a) = K^^^^paPa-<, 


(44) 


where we have defined the following function = fi“(a;“(a))5[a;“(a) — a;“( 6 )], which arises from the definition of V in Eqn. 

by taking the sum over the Bohr frequencies w. The Dirichlet matrix therefore has the same decomposition over cosests 
[ 70 ], syndromes a and single qubit Pauli’s aj. This provides an easy approach to finding bounds on the support numbers tj. 


The dual vectors | a.yg ) are orthogonal if they come from different cosets and thus the full Dirichlet matrix is a direct sum over 
the cosets [ 70 ]. We can therefore write for the matrix 


pQ!D cQ!^ 

S' ^ Vt/ S.70 ’ 
[70] 


(45) 


where each of the matrices can be decomposed into a sum of positive semi-definite two dimensional matrices. We 

can then support the two dimensional matrices from the heat-bath generator by those of the Davies generator, reducing the 
complexity of the problem greatly. 

The matrices and can be decomposed as follows 

Sv7n ~ 2 Qj. 70 '^i,70 ” 2 

a ctj a OLj 

where the matrices Qj^jg{a, aj) and aj) are now two dimensional matrices in the space defined by | and | ). 

To simplify notation we also define the vector I — |a7o) ^a.7|a7(J)- 

We first consider the matrices Qj^-yg{a, aj): It can be verified easily that the functions Gj{a) satisfy the identity Gj{a°‘3^ = 
Gj{a)e^^°‘Ba) for all at site j. This implies that G 3 (a“-> )Pa = Gj{a)pa‘^j so that (a) = B^.^{a°‘). A simple regrouping 
of the sum in Eqn. (l42l i yields Eqn. (l46l l when we define 









11 


Qj,^{a,a) = i y/paPa->\a.y){a-^\ 

+ g “ Gj{a?))^ V^'a“Pa“T'|a“)(a“| 

+ iG,(a“)G,(a“'^)y^;^|-“J(-“J. (47) 

The matrices Dj^jg{a, aj) for the Davies generator can be brought into a similar form: Since the rate function satisfies the 
KMS-condition: i.e. /i“(—w) = h°‘{uj)e~^‘^, we have that h°‘{uj°‘{a))pa = /i“(w“(a“))paa so that B^^{a) = i3^^(a“). An 
important difference with Eqn. ( |47] | is that we take an average over aj = 0 ... 3 for all diagonals and distribute these weights 
equally among the Dj^jg{a, aj). We will explain why this is necessary later. This leads to matrices 

1 ^ 

Djn{a,a) = g X! (^“(w“(a)) + /i“(w“(a^))) (1 - S[uj°‘{a) - uj°‘{a'^)]),/p^\aj){a^\ 

a—0 
1 ^ 

a—0 

+ /r“(w“(a))(5[a;“(a) - w“(a'^)]VPaPa7|-“^)(-“J. (48) 

Note, that the 5 function still appears in this definition. It is therefore important, that we perform a case analysis when comparing 
the matrices. 


We are now finally in the position to state a simple comparison between these two matrices. We find that we can always bound 


max 

a 


G^(a) . 

I— Dj^-f{a, ot) > Qj^j{a,cx), 


(49) 


for all j, a, j and a, so that lifting this constant out of the sum, implies the theorem. To this end we need to consider three 
different cases. We group them according to the values of oj°‘{a). 

1. The first case we consider is when w“(a) = Lo°‘{a^) for all a = 0 ... 3. In this case we have that — a;“(a^)] = 1 

as well as Gj{a°‘) = Gj{a°‘'^) for all a. Therefore, Dj^^{a, a) and Qj^j{a, a) simplify greatly since the purely diagonal 
contributions vanish and we have that 


Dj^^{a,a) = h°'{uj°‘{a))y/paPa'r\-c 


,)K 


and 




)(-“ 

,, / V a-Y 


(50) 


Since the resulting matrices are up to proportionality the same projectors, we only need to compare the coefficients to 
obtain the estimate Tj > maxa^a ^ — ma-Xa 2. The last inequality is obtained by minimizing the 

numerator and the denominator independently. 


2. Now let us assume that all uj°‘{a) ^ so that none coincide. From this we have that 5[a;“(a) — a;“(a^)] = 0 and 

Of) is purely diagonal. Note that the diagonals can be bounded by 4 X]a=o (^“(w“(o)) + /i“(u;“(a'''))) y/paPa~< > 
^min'\/ PaPa'^ and g {G j{jX ) Gj{(l ^)) PaPa^ — g aiclXj ^ G j PaPa^ i SO that 


Cr) ^ ^min'x/ PaPa'^ Ila,a“ 2 ^min'x/ PaPa'^ ^lla,a“ 2 ^ a^}( a II and 

Qj,^ia,a) < imaxG|(a)VPaPaTla.a° + )(-“ 

o 4 T T 


(51) 


as semidefinite inequalities. Since we have brought the matrices in a similar form again, we only need to read off the 
coefficients to obtain the two estimates r, > max„ for the first term and r, > max„ for the second. Both 

J — “ 8/lmin J — “ tlmln 


bounds are again consistent with Tj > maxa 2 


G(a) 


3. The final case to consider is when the Bohr frequencies u!°‘{a),uj°‘{a^) coincide for some a while for others they do not. 
For the a, a where the Bohr frequencies do not coincide, the analysis is similar to case 2 studied above (up to a factor of 
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2 which yields the lower bound). For the a, a with = a;“(a''') we have a different splitting, where 


(tr, rr) ^ ^^minv/ PaPa"^ lla,a“ 2 ^minv/ PaPa'^ \ a.y ) ( ttnd 

Qj,jia,a) < ^m£ixG'^j{a)y/PaPa-rla,a‘^ + ^Gi(a“)G'i(a“^)VPaPa7|-“J(-“^ 


(52) 


It now becomes clear why the redistribution and averaging over a in the definition of the diagonals of Dj^^{a, a) has 
become necessary. Had we not performed the sum, the diagonal contribution in Dj^^{a,a) would have disappeared 
and the range of a) could not be supported by Dj^^{a, a) making it impossible to find a finite support number. 

However, with this splitting we also find that Tj > maxa 2 ^^-^ suffices. 

^min 

From the analysis of these three different cases we have that the bound in Eqn. (|4^ holds. If we now define r = 

maxj;a 2^f^, we have that for all cosets [ 70 ] and all j in the decomposition, that which then translates to 

the final bound in the lemma. □ 


Bound on k 


We now proceed to estimate the for the heat-bath map under the assumption that the Hamiltonian is of commuting Pauli 
form. Before we estimate the constant k, we state bounds on the in terms of suitable norms of 7 ^ and dk"ij which will then 
be bounded in a later step. 

We denote by Nj the set of qubits that share a stabilizer generator gk with site j. That is I G Nj if for some gk both 
l,j € supp(pfc). Moreover, we denote by Bi the ball that is generated by a chosen order. That is, we order all sites I G Nj, with 
the constraint that j itself is the smallest element. Then Bi denotes all elements that are smaller or equal to I with respect to this 
ordering. 

Lemma 4 Let Q denote the heat Bath generator of Eqn. then we have: 

\\[dk,Cj]{f)\\<Y.<M(f)\\^ ( 53 ) 

i 


where 



^WljW '^seB{j)f I|(^s7jll 
0 


if k G Nj 
otherwise. 


(54) 


Note that B{j)i refers to the labels I that are in the chosen ordering in B(j)i C Nj smaller than s < 1. The sum is taken over 
the complement B{j)i = A\B{j)i. 


Proof. It follows immediately from the definition of Qj = Ej [/] — 5 in terms of the non-commutative expectation value that 
\dk,Cj]{f) = 0 whenever k ^ Nj. This directly implies the second line in the definition of the ail- To obtain the bound for the 
case where k G Nj we proceed by recalling the definition of the heat-bath generator 

e<t) i E (Nf^i - ■ (55) 

j,a ^ '' 

where each of the Lindblad operator is defined as 7l“ = 7 j,(t“ with 7 ^ = ^{a)jP{a). We only need to consider the local 

summands in the expressions for Every Lindblad generator can conveniently be written as the sum of two commutators 

= 5 E7(/) = \ E + \Ny-N) ■ (5<') 

a. a. ^ ^ 

Erom this we can directly evaluate the commutator for every summand 

/:“](/) = ^ {[AfMk (/)]7l“ + AfitTk (/) ,71“]) - i (trfc {[Afj]Af) + tr^ (^if [/,7l“])) 


(57) 
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We now proceed by applying a trick in Ref. to evaluate the commutator. Note that [g — tr (g), f — tr (/)] = [g, f] for any 
observables /, g. Then, we can choose an ordering of lattice sites starting at site i, i.e. i is the smallest element in to reach 
the following operator identities, 

N N-1 

g-tT{g)=Y^ {di^g) and / - tr (/) = ^ , (58) 

m—1 n—0 

where we denote the partial trace over the sites by (•) and the complement by [/i, -. -,= 

Note that the shelling of the space works in the opposite direction for / and g so that 

N N-1 

m—1 n—0 
N-1 N 

= Z Z ('9u 5) : ('9i„+i/)] • 

n—0 m>n+l 

(59) 


In the last line we used the fact that partial trace yields an identity on the subsystem so the commutator 

(i9i„+i/)] = 0 for all m < n + 1. This commutator indetity can now be used to bring the 
terms in Eqn. (fSTl i into a suitable form and apply norm bounds from there. We apply Eqn. ( |59] | to Eqn. (fSTl) and obtain 


N-1 N 


[4,7“I(/)=E E 





n=0 m>n+l 









- (trfe 


..../m]' {^ImAf'j ,tr[Z^ 


A° 


+trfe (^Af ). 


(60) 


The operator norm for two matrices is both sub multiplicative, (||(?/|| < ||5||||/||), and contractive under the partial trace 
(||trs (g) II < II 5 II). Eurthermore the triangle inequality can be used to bound || [/, p]|| < 2||/(;||. Hence, for I ^ j 


\\diAf\\ = \\diAf\\<\\daN (61) 

since 9i(7jO'“) = (dijj)a^ because Pauli operators acting on different sites commute. Moreover we have that ||cr“|| = 1. The 
norm bound for djA'j is not so direct, since the action of the Pauli matrix cr“ and dj do not commute any longer. However, we 
observe that this occurs only for the first summands in m which are of the form trjjjo {djA^^ = tr[j]c (^“) — tr (^^)- Since 
this term appears in the commutator and tr (^“) oc 1, we have that this term vanishes in the commutator. Therefore we can 
immediately state the bound 


N-1 N 

Z (62) 

n—0 m>n+l 


Now recall that II[5fe,£j](/)|| < i Set II '('j ](/)ll to the triangle inequality. Relabeling the summands s = Im and 
I = In+i, we are left with the estimate 


which yields the desired bound. 


N / N 


||[5fc,Q,](/)||<Z Z4|l7,IIP.7,ll mi 


1=1 \ S>1 


(63) 


□ 


Since we now have an estimate for the aj.; - constants in terms of the matrices, we only need to evaluate the operator norms 
of and dmlj. In principle the sum in the definition of the has to be taken over the full lattice with the exception of a ball 
Bi that surrounds the site j. However, it will become evident that the summands o-ki will vanish, whenever I ^ Nj any longer. 
This will follow from the property that \\dm'lj || = 0, once m ^ Nj, since the sum is empty when the ball Bi is larger than Nj. 
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Proposition 5 The condition number k is bounded by 


3(|iV,|-l) 
nip) = max --=- - - 

4 


Oii—O iGNj m>l 


Y. - 1 ), 


(64) 


where 


e{m,j) = max Y Jkek{aj)ek{Tm)- (65) 

k 

Remark: Before we proceed with the proof, let us note that a simpler expression can be obtained, which does not depend on 
the syndromes e{a). The only relevant model dependent parameters are the number of qubits | Nj \ that are connected to the qubit 
j through adjacent stabilizers, and the number of stabilizers that are shared between adjacent qubits m, j denoted by fl S'™!- 
We have that. 


k(/ 3) = max3(|iV,-| - ^ ^(g2/3J|S,-nS,| _ 

^ l£Nj k>l 

which follows from the upper bound on e(m, j ) in terms of the shared number of stabilizers. Note that we also gave a bound on 
1 Set =0 by the smallest summand. The smallest summand can simply be bounded by the exponential of the total 

number of stabilizers \Sj \ that are supported on site j. 


Proof. The central ingredient for obtaining bounds on the stems from finding suitable norm bounds on the pj and dmlj- 
Note that the operator norm of the matrix pj is evaluated immediately, since the projectors P{a) are mutually orthogonal, so 
that 


IITjII =maxGj(a) 



E 


^-2pj\e{aj)\ 


- 1/2 


(67) 


Moreover, we can estimate the norm of dipj by the following sequence of steps 

3 


119m7/II = max 


- ^(G,(a)-G,(a-)) 


rm=0 


< -max|Gj(a) - Gj(a^”*)| 


( 68 ) 


where we have made use of the fact that we can bound every summand by the maximal except for r^, = 0, which always 
vanishes. Now it is easy to see that 


max|G,(a) — Gy(a'^'")| < || 7 ,jl max 

a,Tm CL-,'Tm 




(69) 


We now assume that ) — 1 so we can drop the absolute values. If this is not the case we can map a —?► o’” and flip the 

sign so that the assumption holds again. This is possible since we are considering the maximum over all a. So the two cases are 
in fact symmetric. Bounding both terms for every summand in Q we can see that we always have 


max 

a,Tm 


G,(a--) 


~ TmyOij 


(70) 


since the exponential is a monotone function the bound is nothing but exp(2/3e(m, j)). Hence we obtain the upper bound to the 
norm 


||9m7ill< (71) 

from this bound it can be seen easily that ||9i7jjl = 0 ’ whenever I ^ Nj as would be expected. This holds, since e{l,j) = 0, 
when the sites l,j do not share a stabilizer. Hence using the formula from lemma|4]for the we have 

r 3||7if EmeivAB, - 1) : if fc G TV, 


0 


otherwise. 


( 72 ) 
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Now let us recall that the constant k(/ 3) is determined in therms of the sum 

«.(/3)<max^^<,=max(|iV,|-l)^3||7,f ^ - l) . (73) 

^ kejc I ^ l^Nj mGNj\B, 

The sum in k vanishes for summands k ^ Nj as already observed in lemma|4] Since the summands do not depend explicitly on 
k, we just estimate this sum in terms of | | — 1 because j, the center, is excluded. Moreover, we have observed that dmlj = 0, 

whenever m ^ Nj. This leads to the fact the vanish for all I ^ Nj because the set Nj\Bi is empty. This justifies the 
restriction of the summands in the second equality. Now recall that the norm || 7 jH is given by Eqn. dbTl i so that we are left with 
the expression as stated in the lemma. □ 


Gap from ergodicity 


In this section we show that strong ergodicity of the semigroup (Eqn. (fTTT il implies that the generator is gapped. Indeed, 
Lemma 6 Let C be the generator of a primitive reversible quantum dynamical semigroup. If there exist constants c, k such that 

l|e*^(/)-tr(p/)|U<c|||/|||e-*(i-'=), (74) 

for all observables f, then the generator has a gap which is lower bounded by 1 — k. 


Lg norms. A non-commutative Cp norm is defined 


Proof The proof hinges on the use of interpolation theorems for Lp 
with respect to a full rank density matrix p > 0 as 11/| ^ = tr (see Refs. for more details). The 

relevant inner product is (/, g)p = tr ■ We will first show that 


e*^-E°°||i_i,p= ||e*^-E’ 


tc 


(75) 


where E°° = limt_).oo Indeed, from the variational characterization of norms, and reversibility of the semigroup. 


i.p= sup ||(e‘^- 

-E°°)(/)l|i.P 


ll/lli 

,P<1 



= 

sup 

(p, (e*^ - E“)(/)) 

llslU 




= 

sup 

(/,(e*^-E-)(p)), 

llsllo. 

= <1HI/Ili,p<i 


= sup ||(e‘^- 

E“)(5)IIoo 


Hallo. 

= <1 




=||e‘^-E“||oo-oo 


where we have used that E°“ and are reversible with respect to the inner product (/, g)p = tr ^p 2 pp^ . Now, invoking 
an extension of the Riez-Thorin theorem CSll (see Lemma|7]i, we get for any linear operator T, 


||T(/)||2.p<2y/2iViiV„ 


2,P5 


(76) 


whenever ||T||i_i,p < Ni and ||r||oo-oo < N^o- 

Therefore, setting T := — E°°, and observing that 

\PHf) - tr(p/) lU < c|||/|||e-(i-'‘)‘ < 2c|Ay| \\f\\^e-^^-^^^ (77) 

where |A/| is the support of observable / and 

||e‘^(/) - tr {pf) ||i,p < c|A/| (78) 

we get 

l|e‘^(/)-tr(p/)||2,p<2c|Ay|||/||2.pe-(i-'^)*. 


( 79 ) 
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Then, by picking the operator f\ such that f\ — tr {pf\) is the eigenvector corresponding to the spectral gap A, we get 

l|e‘^(/A) -tr(p/A) II 2 .P = I I/a - tr (p/a) (80) 

<c'|A/| ll/Alb.pe-^i-'^)* (81) 

This inequality must holds for all f > 0, hence A > 1 — k. □ 


Lemma 7 (Ref. Il28ll ) Let T be a bounded linear operator, and assume that there exist constants Np, Nq such that for p < q, 
and any (hennitian) operator f we have 


ll"^(/)llp,p ^ ^PWJ iip.p 
then for any p < r < q we get 


and ||T(/)||g,p < 


q,p SI \\q,P 


||T(/)||..p<7iVpAf^^uii..p, 


with 


5 = 


piq-r) 

r{q-p) 


and 


7 = 2 


r{q-p) 

{r-p){q-r) 


(82) 


(83) 


(84) 






